Immune-related 3-lncRNA signature with prognostic connotation in a multi-cancer setting.

BACKGROUND: Advances in our understanding of the tumor microenvironment have radically changed the cancer field, highlighting the emerging need for biomarkers of an active, favorable tumor immune phenotype to aid treatment stratification and clinical prognostication. Numerous immune-related gene signatures have been defined; however, their prognostic value is often limited to one or few cancer types. Moreover, the area of non-coding RNA as biomarkers remains largely unexplored although their number and biological roles are rapidly expanding. METHODS: We developed a multi-step process to identify immune-related long non-coding RNA signatures with prognostic connotation in multiple TCGA solid cancer datasets. RESULTS: Using the breast cancer dataset as a discovery cohort we found 2988 differentially expressed lncRNAs between immune favorable and unfavorable tumors, as defined by the immunologic constant of rejection (ICR) gene signature. Mapping of the lncRNAs to a coding-non-coding network identified 127 proxy protein-coding genes that are enriched in immune-related diseases and functions. Next, we defined two distinct 20-lncRNA prognostic signatures that show a stronger effect on overall survival than the ICR signature in multiple solid cancers. Furthermore, we found a 3 lncRNA signature that demonstrated prognostic significance across 5 solid cancer types with a stronger association with clinical outcome than ICR. Moreover, this 3 lncRNA signature showed additional prognostic significance in uterine corpus endometrial carcinoma and cervical squamous cell carcinoma and endocervical adenocarcinoma as compared to ICR. CONCLUSION: We identified an immune-related 3-lncRNA signature with prognostic connotation in multiple solid cancer types which performed equally well and in some cases better than the 20-gene ICR signature, indicating that it could be used as a minimal informative signature for clinical implementation.

Network-based identification of key master regulators associated with an immune-silent cancer phenotype.

A cancer immune phenotype characterized by an active T-helper 1 (Th1)/cytotoxic response is associated with responsiveness to immunotherapy and favorable prognosis across different tumors. However, in some cancers, such an intratumoral immune activation does not confer protection from progression or relapse. Defining mechanisms associated with immune evasion is imperative to refine stratification algorithms, to guide treatment decisions and to identify candidates for immune-targeted therapy. Molecular alterations governing mechanisms for immune exclusion are still largely unknown. The availability of large genomic datasets offers an opportunity to ascertain key determinants of differential intratumoral immune response. We follow a network-based protocol to identify transcription regulators (TRs) associated with poor immunologic antitumor activity. We use a consensus of four different pipelines consisting of two state-of-the-art gene regulatory network inference techniques, regularized gradient boosting machines and ARACNE to determine TR regulons, and three separate enrichment techniques, including fast gene set enrichment analysis, gene set variation analysis and virtual inference of protein activity by enriched regulon analysis to identify the most important TRs affecting immunologic antitumor activity. These TRs, referred to as master regulators (MRs), are unique to immune-silent and immune-active tumors, respectively. We validated the MRs coherently associated with the immune-silent phenotype across cancers in The Cancer Genome Atlas and a series of additional datasets in the Prediction of Clinical Outcomes from Genomic Profiles repository. A downstream analysis of MRs specific to the immune-silent phenotype resulted in the identification of several enriched candidate pathways, including NOTCH1, TGF-$\beta $, Interleukin-1 and TNF-$\alpha $ signaling pathways. TGFB1I1 emerged as one of the main negative immune modulators preventing the favorable effects of a Th1/cytotoxic response.

A modeling framework for embedding-based predictions for compound-viral protein activity.

MOTIVATION: A global effort is underway to identify compounds for the treatment of COVID-19. Since de novo compound design is an extremely long, time-consuming and expensive process, efforts are underway to discover existing compounds that can be repurposed for COVID-19 and new viral diseases.We propose a machine learning representation framework that uses deep learning induced vector embeddings of compounds and viral proteins as features to predict compound-viral protein activity. The prediction model in-turn uses a consensus framework to rank approved compounds against viral proteins of interest. RESULTS: Our consensus framework achieves a high mean Pearson correlation of 0.916, mean R2 of 0.840 and a low mean squared error of 0.313 for the task of compound-viral protein activity prediction on an independent test set. As a use case, we identify a ranked list of 47 compounds common to three main proteins of SARS-COV-2 virus (PL-PRO, 3CL-PRO and Spike protein) as potential targets including 21 antivirals, 15 anticancer, 5 antibiotics and 6 other investigational human compounds. We perform additional molecular docking simulations to demonstrate that majority of these compounds have low binding energies and thus high binding affinity with the potential to be effective against the SARS-COV-2 virus. AVAILABILITY AND IMPLEMENTATION: All the source code and data is available at: https://github.com/raghvendra5688/Drug-Repurposing and https://dx.doi.org/10.17632/8rrwnbcgmx.3. We also implemented a web-server at: https://machinelearning-protein.qcri.org/index.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Ancestry-associated transcriptomic profiles of breast cancer in patients of African, Arab, and European ancestry.

Breast cancer largely dominates the global cancer burden statistics; however, there are striking disparities in mortality rates across countries. While socioeconomic factors contribute to population-based differences in mortality, they do not fully explain disparity among women of African ancestry (AA) and Arab ancestry (ArA) compared to women of European ancestry (EA). In this study, we sought to identify molecular differences that could provide insight into the biology of ancestry-associated disparities in clinical outcomes. We applied a unique approach that combines the use of curated survival data from The Cancer Genome Atlas (TCGA) Pan-Cancer clinical data resource, improved single-nucleotide polymorphism-based inferred ancestry assignment, and a novel breast cancer subtype classification to interrogate the TCGA and a local Arab breast cancer dataset. We observed an enrichment of BasalMyo tumors in AA patients (38 vs 16.5% in EA, p = 1.30E - 10), associated with a significant worse overall (hazard ratio (HR) = 2.39, p = 0.02) and disease-specific survival (HR = 2.57, p = 0.03). Gene set enrichment analysis of BasalMyo AA and EA samples revealed differences in the abundance of T-regulatory and T-helper type 2 cells, and enrichment of cancer-related pathways with prognostic implications (AA: PI3K-Akt-mTOR and ErbB signaling; EA: EGF, estrogen-dependent and DNA repair signaling). Strikingly, AMPK signaling was associated with opposing prognostic connotation (AA: 10-year HR = 2.79, EA: 10-year HR = 0.34). Analysis of ArA patients suggests enrichment of BasalMyo tumors with a trend for differential enrichment of T-regulatory cells and AMPK signaling. Together, our findings suggest that the disparity in the clinical outcome of AA breast cancer patients is likely related to differences in cancer-related and microenvironmental features.